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ABSTRACT. 

The Sun’s gravitational field deflects the apparent positions of close objects in accordance with the 
formulae of general relativity. Optical astrometry is used to test the prediction, but only with the stars 
close to the Sun and only during total Solar eclipses. Geodetic Very Long Baseline Interferometry (VLBI) 
is capable of measuring the deflection of the light from distant radio sources anytime and across the whole 
sky. We show that the effect of light deflection is equivalent to the gravitational delay calculated during 
the reduction of VLBI data. All reference radio sources display an annual circular motion with the 
magnitude proportional to their ecliptic latitude. In particular, radio sources near the ecliptic pole draw 
an annual circle with magnitude of 4 mas. This effect could be easily measured with the current precision 
of the geodetic VLBI data. 


1. INTRODUCTION 

Very Long Baseline Interferometry (VLBI) is capable of measuring precise group delays - the difference 
in arrival times of radio waves at two radio telescopes (Schuh & Behrend, 2012) from distant extragalction 
radio sources (quasars). Accurate positions of these radio sources are obtained with an accuracy of 40 
microarcsec (Ma et ah, 2009). 

Gravitational time delay caused by the Solar gravitational field is calculated during the reduction of 
geodetic VLBI data (Shapiro, 1964, 1967). The conventional formula for calculating gravitational delay 
is formulated in terms of the positions of the radio telescopes within the barycentric reference frame of 
the Solar System (Kopeikin, 1990; Eubanks et ah, 1991; Soffel et ah, 1991; Klioner, 1991), rather than 
the baseline length between the radio telescopes. 

We propose an alternate gravitational delay formula using a Taylor series expansion. We show that 
the conventional formula can be split into a sum of several terms, and the major term links the gravi¬ 
tational delay and the well-known formula for the light deflection angle. The light deflection angle can 
be considered equivalent for all baselines and estimated for each radio source at times of interest. We 
develop a new approach to probe the formula for the light deflection angle at an arbitrary elongation 
from the Sun. Finally, estimates of the light deflection for several reference radio sources based on VLBI 
observations in I99I-200I are presented. 

2. GRAVITATIONAL DELAY VS LIGHT DEFLECTION 

The light deflection angle a at an arbitrary elongation Q from the Sun, is given by 

_ (7 -I- l)GM sin0 
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where G is the gravitational constant, M is the mass of a gravitational body, c is the speed of light, r is 
the distance from the Earth to the Sun, 7 is the parameter of the Parametrised Post-Newtonian formalism 
(PPN), 6 is the elongation angle (i.e. the angular distance between the Sun and the gravitational body, 
as measured by the Earth observer) (Shapiro, 1967; Ward, 1970). The conventional gravitational delay 



is calculated as follows 
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The VLBI total delay model also comprises a term due to the transformation from the barycentric to the 
geocentric reference frames. 
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It can be shown that the three formulae are approximately linked as (Titov & Girdiuk, 2015) 


'^GR ~ Rgrav Tcoord — G sin (p COS A (4) 

where angles p, 9, A are linked by the standard spherical triangle formula 

cos^/i = — cos7)cos0 — sinTJsindcos^, (5) 

Converting the three dot products with angles ip, ■0, 9 (Fig 1), where 

(& ■ s) = |6| cos(p, (5 ■ rj) = |5||r5| COS0, (rj • ^ = -|r5| cos0 (6) 



Figure 1: Angles p, ip, 9 and A, originated by posi¬ 
tion of gravitational mass (B), quasar (Q) and base¬ 
line vector (b). If the Sun plays the role of the grav¬ 
itational mass, then, the point B on Fig 1 is also 
the position of the Solar system barycentre. 


More accurately, the total contribution of the post-Newtonian (PN) general relativity effect to the 
total VLBI delay is given by (we recall now 7=1 for the case of general relativity) 
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The two additional terms appear because the elongation angles 9i and 02 between the direction to the 
observed radio source and the centre of the gravitational body, as measure from each radio telescope, 
are not equivalent. The parallactic effect for a baseline of 6,000 km is about 8 arcsec, and the deflection 
angles ai and q ;2 for the ’’first” and the ’’second” radio telescope differing. At a large elongation the 
additional terms are negligible, so in the small angle approximation this effect should be considered. 

TGR = c y’(cos A — sin p cos 2A) = q ;2 7 p{cos A — sin p cos 2A) ( 8 ) 

where a 2 is the classical deflection angle for light propagated through a gravitation field (Einstein, 1916) 
for the second station a 2 = and i ?2 = ^"2 cos 02 - In an alternative form 
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where ^ the additional deflection angle 


3. OBSERVATIONAL DATA 

We processed the available VLBI data between 1991 and 2001 using the OCCAM software (Titov et 
ah, 2004). A fraction of well-established reference radio sources were assigned as astrometrically unstable, 
i.e. their coordinates were not fixed. The observational group delays were approximated by the theoretical 
values, and the 0-C (observed minus calculated) differences were used to obtain a set of estimates of the 
standard daily parameters (Earth orientation, station positions, wet troposphere delays and gradients, 
etc.) and corrections to the selected radio source coordinates (right ascension and declination) regardless 
of its elongation from the Sun. The consensus gravitational delay (3) and the coordinate term (4) were 
not applied for the calibration of the group delay for the radio sources. The ionosphere fluctuations were 
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Figure 2: The variations of the coordinates of the radio sources 1606+106 (left) and 0229+131 (right) in 
1991-2001. A smaller minimum angle 6 for 0229+131 results in larger size of the circle. 


calibrated in the conventional way, and, as it will be shown later, this makes possible VLBl observations 
of radio sources in the range of 1°.5 to 3° from the Sun. 

We used the Least Squares collocation method to estimate the wet troposphere fluctuations with the 
OCCAM software. The mutual correlation between observables are introduced in the data adjustment 
process. The difference between the VLBl estimates of the wet troposphere delays and independent 
radiometer data is typically within 3-6 mm, or 10-20 ps (Titov & Stanford, 2013). Thus, the impact of 
the wet troposphere delay on the astrometric light deflection angle estimated near the Sun is negligible. 
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Figure 3: Variations of the daily estimates of the PPN parameter 7 for radio source 0229+131 in 1991-2001 
at the different scale of angle 6. 


Fig 2 shows the variations of the coordinates of two selected radio source over 10 years (1991-2001). 
Each point on the plot corresponds to a daily position estimated from a single 24-hour session. About 
10 annual circles are drawn by each radio source on the plots. The catalogue (non-deflected) positions 
are at the reference origin, and the total deflection angle is equal to distance to the point from the 
reference origin. The ecliptic latitude of the four sources are different, therefore, elongation angles vary 
over different ranges. As a result, the deflection angle does not change essentially for the former radio 
sources (the reference origin is in the middle of the ring) with a magnitude about 4 mas. In contrast, for 
the latter case, the deflection angle varies from zero to 0".2 (for the minimum elongation of 3°), and the 
reference origin lies near the edge of the ring. 

The variations of the deflection angle are easily converted with formula (1) to the variations of the PPN 
parameter 7 . Fig 3 shows the variations of the deflection angle as a function of 6 for the reference radio 
source 0229+131 at different scales. This radio source is almost on the ecliptic plane (ecliptic latitude 
e ~ -1°.5), and the elongation angle is limited by the technical capability of observing near the Sun. 
On 29 April, 1997 (session code 97APR29XE) the mean elongation angle was 1°.55, and the measured 
deflection reached a maximum among all selected radio sources, a = 0".3100 +/- 0".0009. However, this 







estimate was obtained from only 7 observations. The best accuracy was achieved on 1 May, 1996 (session 
code 96MAY01XO) with 122 observations. The deflection angle estimate is a = 0".1588 +/- 0".0002, 
and the relative accuracy is 0.0012, in spite of a larger mean elongation angle for this date. 

4. CONCLUSION 

The effects of general relativity are explicitly contained in both components of the total VLBI delay 
model - the gravitational and geometric delays. While the former component uses the individual barycen- 
tre positions of the radio telescopes to calculate the effect, the latter component is expressed in terms of 
the baseline between the radio telescopes. Coupling between both parts has not been investigated until 
now. 

The Shapiro effect as measured by radars, and the deflection of light measured with traditional 
astronomical instruments are considered two independent tests of general relativity. In this paper we 
show that the total group delay model joins the two tests within one observational technique - geodetic 
VLBI. The gravitational delay that traditionally originated from the Shapiro effect is linked to the light 
deflection angle. Therefore, the two approaches, VLBI delay and angular, are absolutely equivalent. 

For almost all realistic situations this angle does not depend on the baseline length, thus, a standard 
geodetic VLBI interferometer acts as a traditional astronomical instrument. In addition, the coordinate 
term explicitly presented in the conventional geometric delay model ceases to exist because it is compen¬ 
sated by the same effect in the gravitational delay with the opposite sign. Thus, the proposed alternative 
version of the general relativity contribution to the total VLBI group delay model is free from coordinate 
effects. Therefore, the two approaches, time delay and angular deflection, are absolutely equivalent. 

The final equation of the general relativity contribution also comprises two smaller terms which are 
significant at very small angular separation between the deflecting body and distant radio source. These 
terms may be considered as an increment in the light deflection angle due to the additional time delay 
for propagation of the light from station I to station 2. This effect become significant at 4 > O.I. 
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